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Abstract: A comparative study of the exclusive rare B c — > D*£ + £~ (£ = /i,r) decays 
has been made in the minimal supersymmetric models (MSSM) and the SUSY SO (10) 
GUT models. In this context, various physical observables such as branching ratios (BTZ), 
forward-backward asymmetries (Afb), lepton polarization asymmetries (Pl,n,t) and he- 
licity fractions (/l,t) of D* meson by using the the QCD sum rules form factors have 
been investigated. It is found that the SUSY effects are characteristically prominent to 
that of the SM values for these observables. For instance, in SUSY I and SUSY II, the 
forward-backward asymmetry does not cross zero which is mainly due to the same sign of 
the C e 7 n and Cg ff Wilson coefficients. Similarly in SUSY SO (10) GUT models due to the 
complex nature of the new Wilson coefficients - corresponding to the new operators arising 
due to the contribution of neutral Higgs bosons (NHBs) - the above mentioned observables 
are sizably affected. Therefore the analysis of said observables in charmed semileptonic B 
meson decays can put some stringent constraints on the parameter space of SUSY variants 
and can serve as a windowpane to look beyond the SM. 
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1 Introduction 

The standard model (SM) of particle physics is considered to be one of the most successful 
theory of the twentieth century but we can not accept it as the ultimate theory of nature 
since there are many open questions beyond the scope of the SM to be addressed. These 
questions include: gauge and fermion mass hierarchy, matter-antimatter asymmetry, num- 
ber of generations of quarks and leptons, the nature of the dark matter and the unification 
of fundamental forces. Therefore, we need to search new physics (NP) beyond the SM that 
may help us to answer these open problems of the SM. One of the plausible extension of 
the SM is supersymmetry (SUSY) [1] which is a promising candidate of NP and can shed 
light on some of the issues mentioned above [2]. For example, SUSY solves the hierarchy 
problem of the SM by the cancellation of the quadratic divergences in the radiative loop 
corrections since in SUSY fermions and bosons contribute with opposite sign [3, 4]. It 
is also considered that the lightest supersymmetric particles (LSP) are the dark matter 
candidates because they are stable and interact very weakly with the ordinary matter. It 
is also an important ingredient in the superstring theory which is a suitable candidate for 
the unification of all the known forces including gravity. 

There are two ways to search for the SUSY, one is to discover the SUSY particles 
(sparticles) at high energy colliders directly [5] and the other is to search for its effects 
through indirect methods. Since no sparticle has been seen so far, they must have higher 
masses than their SM partners, implying that the SUSY is a broken symmetry. The other 
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option is to investigate the SUSY through the indirect searches where SUSY can show 
its footprints. In this regard, the processes involving the flavor changing neutral currents 
(FCNC) are the elegant way to probe the indirect searches of the SUSY. So the rare decays, 
involving FCNC, induced by b — > s(d) transitions provide a windowpane to look for the 
physics beyond the SM. This lies in the fact that FCNC processes are not allowed at the 
tree level in the SM but can occur at the loop level through Glashow-Iliopoulos-Maiani 
(GIM) mechanism [6] thereby making them quite sensitive to the indirect searches of the 
SUSY [7, 8]. The sparticles can contribute to FCNC transitions through the quantum loop 
due to R-parity conservation, hence, making these processes a handy tool for studying the 
possible effects of sparticles in the rare decays [9—11]. This gives us a good motivation 
to study the rare decays both theoretically and experimentally to look beyond the SM 
especially in the LHC era [12, 13]. 

Since FCNC processes are only possible at the loop level and have small branching 
ratios as compared to the tree level transitions but they have implications for the results 
obtained in the ongoing experiments after the observations of the rare radiative b — > s'j 
decay at CLEO [14]. Since then there have been many studies on the rare semileptonic, 
radiative and leptonic decays of B u ^ s mesons induced by the FCNC transitions of b — > s(d) 
[15] . These studies will be even more complete if one considers similar decays of the charmed 
B mesons (B c ). The charmed B c meson is a bound state of two heavy quarks, the bottom b 
and the charm c, and was first observed in 1998 at the Tevatron in Fermilab [16] . Because 
of the two heavy quarks, the B c mesons are rich in phenomenology compared to the other B 
mesons. At Large Hadron Collider (LHC) the expected number of events for the production 
of B c mesons are about 10 8 — 10 10 per year [17-22] which is a reasonable number to work 
on the phenomenology of the B c meson. In the literature, some of the possible radiative 
and semileptonic exclusive decays of B c mesons such as B c — > (p, K*,D*, B*) 7, B c — > £v^, 
B c -> B*£+£-, B c -> D\tv, B c -> D* sQ £ + £~ and B c -> D* s /+£~ have been studied using 
the frame work of relativistic constituent quark model, QCD Sum Rules, the Light Cone 
Sum Rules and the Ward identities [23-33]. 

Theoretically what makes the B c — > D*£ + £~ more important compared to the other 
B meson decays, such as B° — > (K* , K±, p,ir)£ + £~ , is the fact that it can have two types 
of contributions, one is through the penguin (loop level) and the other is due to weak 
annihilation (WA). In the ordinary B meson decays the WA contributions are very small 
and can be ignored as compared to the corresponding penguin contributions. However, for 
the B c meson these WA contributions are proportional to the CKM matrix elements V c bV* s 
and hence can not be ignored as compared to penguin contributions [33]. In perspective 
to the new physics (NP) effects in the presence of weak annihilation (WA) contribution, 
it is shown [34] that for B° — > K*p + p~ the helicity fractions are mildly effected under 
the implication of SM4 while in the case of B c — > D*p + p~ the SM4 effects are very op- 
timistic [35]. Similarly for this observables when tauns are the final state the effects are 
comparable i.e. for B c — > D*p + fi~ the shift in the maximum(minimum) values in the lon- 
gitudinal(transverse) helicity fractions is 0.23 and in the case of B° — > K*t + t~ this shift is 
0.20. With this motivation it is also interesting to study the SUSY effects for the lepton po- 
larization asymmetries in this decay channel which is calculated in the present manuscript. 
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The present study shows that the lepton polarization asymmetries for B c — > D*l + l~ more 
influenced due to the supersymmetric models in comparison of B° — > K*l + l~ [36] and one 
can clearly see from the graphs that in the transverse lepton polarization for muons and 
tauns the NP effects are more prominent for B c — > D*l + l~ . This feature arises due to the 
interference between the WA and the Supersymmetric contribution which is encapsulated 
in the ft and /g auxilary functions. In the study of the exclusive .B-meson decays the 
form factors which are the non-perturbative quantities and are the scalar functions of the 
square of momentum transfer are important ingredients. In literature the form factors for 
B c — > D*£ + £~ decay were calculated using different approaches including the light front 
constituent quark models, the relativistic quark models, the Ward identities and the QCD 
sum rules [25-33, 37, 38]. In this study we borrow the form factors calculated by the QCD 
sum rules approach for our numerical calculations from Azizi et al. [38]. However, in the 
context of hadronic uncertinities which enter through the relevent form factors, it needs to 
be stressed that various asymmetries such as forward-backward asymmetries, polarization 
asymmetries of the final state meson and leptons (which are calculated in the manuscript 
submitted) have almost negligible influence. Therefore, any deviation from the SM in these 
observables is the clear indication of NP. 

It is important to emphasize here that the NP effects manifest themselves in the rare B 
decays in two different ways, one is due to the new contribution to the Wilson coefficients 
and other is due to the appearance of new operators in the effective Hamiltonian, which are 
absent in the SM. The manifestations of the NP due to the SUSY is unique in a sense that 
it modifies the Wilson coefficients as well as it introduces the new operators in the effective 
Hamiltonian. In the present study, the NP effects are analyzed by studying the branching 
ratio BIZ, the forward-backward asymmetry Afb, the lepton polarization asymmetries and 
the helicity fractions of D* meson for B c — > D*£ + l~ decays both in the SM and in different 
MSSM scenarios with special emphasis on the effects of neutral Higgs bosons (NHBs) [39- 
44]. A sizeable deviations to the SM results due to the SUSY effects are observed in the 
above mentioned observables for many rare decays [45-50]. 

It is worth mention that previously an anomaly has been observed by Belle [51] in the 
lepton forward-backward asymmetry (Afb) in the exclusive decay B — > K*n + /j,~ [52-55]. 
It was noticed that in the large q 2 region q 2 > 14 GeV 2 , the Afb(q 2 ) measurements tend 
to be larger than that of the SM predictions, although the behavior was same for the both 
cases. The discrepancy was more severe at the low q 2 region where the SM prediction 
is negative [56], whereas the data favor positive values and Belle claimed that this is the 
clear indication of NP [59] . Contrary to this the LHCb recently announced the first result 
on the lepton forward-backward asymmetry (Afb) in the exclusive decay B — > K* fj> + [i~ 
which is close to the SM predictions. In particular, the position of the zero crossing of 
the Afb(q 2 ) he on the same position as predicted by the SM [9, 57, 58]. This situation is 
really exciting and tighten the screws on scenarios those have not zero-crossing point of the 
Afb(q 2 )- In this context, SUSY which previously become a prime candidate to explain the 
Belle anomaly because in some of the scenarios of MSSM such as SUSY-I and SUSY-II, the 
Afb(q 2 ) do not cross the zero position [45-50] now seems to be ruled out after the LHCb 
data. However, before to say any final remarks about SUSY-I nd SUSY-II, it needs to 
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be required some confirmation through complementary information as LHCb collaboration 
said that they will continue to collect more data and try to see deviation in the observable 
data if there is any NP exist. With this motivation we have predicted the MSSM effects 
in the B c — > D*l + l~ , where I = /x,t, decay not only for the Afb{<1 2 ) but also for the other 
asymmetries which depends on the polarization of the final state particles such as lepton 
polarization asymmetries and the helicity fractions. 

In this paper, we investigate the decay processes B c — > D*£ + £~ in the context of the 
SUSY SO(10) GUT models [60]. Since the effects of the counterparts of usual chromo- 
magnetic and electro-magnetic dipole moment operators as well as semileptonic operators 
with opposite chirality are suppressed by msjm^ and as such become insensitive in the SM. 
However, in the SUSY SO(10) GUTs their effects can be significantly large, since complex 
flavor non-diagonal down-type squark mass matrix elements of 2nd and 3rd generations in 
the RR sector can be as large as 0.5 [60, 61]. Furthermore, 5^3 can m duce new 

operators, the counterparts of usual scalar operators (Qi,2) in the SUSY models, due to 
the NHBs penguins with gluino-down type squark propagated in the loop. The values of 
the relevant new Wilson coefficients for the MSSM and the SUSY SO(10) GUT models are 
collected from Refs. [45, 46]. The NHBs could contribute largely to the inclusive processes 
B — > X s £ + £~ , because part of the SUSY contributions is proportional to the tan 3 /? [62, 63]. 
Subsequently, the physical observables, such as the branching ratio BIZ, forward-backward 
asymmetries Afb, lepton polarization asymmetries and helicity fractions of the final state 
meson, in the large tan f3 region of parameter space in the SUSY models can be quite 
different from that in the SM. Hence, the measurement of these observables will give us 
some hints of the SUSY effects in these decays. 

The paper is structured as follows. In Sec. 2 we present the theoretical framework 
for the decay B c — > D*£ + £~ necessary for the study of the different SUSY variants. In 
Sec. 3 we present the basic formulas for physical observables such as decay rate, forward- 
backward asymmetries Afb, lepton polarization asymmetries and helicity fractions of D* 
meson. The numerical analysis and discussion on these observables is given in Sec. 4. 
Section 5 gives the summary of the results of our study. 

2 Theoretical Framework 

In this section we give the effective Hamiltonian for B c — > D*£ + £~ decays. We can split the 
contributions of the total decay amplitude into two, one is the penguin (FCNC trasitions) 
and the second is weak annihilation contribution. 

2.1 Penguin Amplitude 

At quark level, the semileptonic decay B c — > D*£ + £~ is governed by the transition b — > 
s£ + £~ for which the general effective Hamiltonian in the SUSY SO (10) GUT model, can 
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be written, after integrating out the heavy degrees of freedom in the full theory, as [46] 



H eff ~ 7K V tb V ts 



ac r * 10 



V2 



i=l i=3 

8 i 

^{cMQiM + cW&M} , (2-i) 

i=l 



where Oi(/u) (i = 1, 2, • • • , 10) are the four quark operators and Ci(fi) are the corresponding 
Wilson coefficients at the energy scale fx [39-44] which is usually taken to be the 6-quark 
mass (mb)- The theoretical uncertainties related to the renormalization scale can be re- 
duced when the next to leading logarithm corrections are included. The new operators 
Qi(i = 1,2, ••• ,8) come from the NHBs exchange diagrams, whose manifest forms and 
corresponding Wilson coefficients can be found in [64-68]. The primed operators are the 
counterparts of the unprimed operators, which can be obtained by flipping the chiralities in 
the corresponding unprimed operators. It is worth mentioning that these primed operators 
will appear only in the SUSY SO(10) GUT model and are absent in the SM and the MSSM 
[45]. 

The explicit forms of the operators responsible for the decay B~ — > D*~£ + £~ , in the 
SM and the SUSY models, are 

(2.2a) 
(2.2b) 
(2.2c) 
(2.2d) 
(2.2e) 
(2.2f) 



o 7 = 


e 

16vr 2 


m b (sa^Rb)F^ 


0> = 


e 2 
16vr 2 


m b (sa, lu Lb)F^ 


o 9 = 


e 2 
16vr 2 


{s ltl Lb)I^l 


o' 9 = 


e 2 
16vr 2 




Oxo = 


e 2 
16vr 2 


{s ltl Lb)h^H 


o' w = 


e 2 
16vr 2 


(s 7 „Rb)e^ 7 5 £ 



e 2 



Qi = j^{sRb)H (2.2g) 
Q'i = ^{sLb)U (2.2h) 



e 



2 



Q 2 = ^{sRb)h b l (2.2i) 
Q' 2 = -^2^)h^ (2-2j) 



with L,R= \ (l=F7 5 )- 
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Using the effective Hamiltonian given in Eq.(2.1) the free quark amplitude for b 
s£ + £~ can be written as 



M P ™ G (b ^ s£+r) = -^ VtbV * 

\/2tt 



c e 9 ff ( M ) {s^mhn) + c' 9 eff (is) {s^Rb)(hn) 

+ C w (s^Lb)(£^ 7 5 £) + C[ (s^Rb)(£ 7 ^ 5 £) 

- 2c e 7 ff (p) "^(sio^Rb)hn - 2c' 7 eff (fj.) ^(sia^Lby^n 

q z q z 
+ C Q1 (sRb) (££) + C Q1 {sLb) (££) + C Q2 (sRb) (l^H) + C' Q2 (sLb) (Z 7 5 £) 

(2.3) 

where q is the momentum transfer. Note that the operator Oiq given in Eq.(2.2e) can 
not be induced by the insertion of four quark operators because of the absence of Z-boson 
in the effective theory. Therefore, the Wilson coefficient C\q does not renormalize under 
QCD corrections and is independent of the energy scale fj,. Additionally the above quark 
level decay amplitude can get contributions from the matrix element of four quark opera- 
tors, Ya=i {£ + £~s \ Oi\b) , which are usually absorbed into the effective Wilson coefficient 
Cg (fx) and can be written as [69-75] 

C 9 e// (/i) = CM + Y SD (z, s') + Y LD (z, s'). 

where z = m c /mb and s' = q 2 /m 2 . Ysd{z,s') describes the short distance contributions 
from four-quark operators far away from the cc resonance regions, and this can be calcu- 
lated reliably in the perturbative theory. However the long distance contribution Yld(z, s') 
cannot be calculated by using the first principles of QCD, so they are usually parametrized 
in the form of a phenomenological Breit-Wigner formula making use of the vacuum satura- 
tion approximation and quark hadron duality. The expressions for the short-distance and 
the long-distance contributions Ysd(z,s') is given as 

Y SD (z,s') = h(z,s') [3C 1 (fi) + C 2 (fj,) + 3C 3 (pL) + C 4 ^) + 3C 5 ( ( i) + C 6 (pi)} 
-h(l, s') [4C 3 (aO + 4C 4 ( M ) + 3C 5 (fi) + C 6 (m)] 

-^(0, s') [C 3 ( M ) + 3C 4 (/i)] + 2 - [3C 3 (m) + C 4 (m) + 3C 5 ( M ) + CM] (2-4) 
Y LD (z, s') = J-(3Ci( M ) + C 2 (/i) + 3C 3 (fi) + C 4 (m) + 3C 5 (fi) + C 6 ( M )) 



a 2 



, 9x, 7rr(j -> l + l~)M j . , 

j=i>,i>' j J j 



with 



8 8 4 2 [in 

h(z,s') = — Inz + — + -x (2 + x)|l - x\ 1/2 



iir for x = Az 2 /s' < 1 



l-x-l 

27 ' 9 9 V ' " 1 I 2arctan^= for x = iz 2 /s' > 1 



WO, s = — m his H — i7r . (2.6) 

v ; 27 9 fj, 9 9 v ; 
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Here Mj(T t ° t ) are the masses (widths) of the intermediate resonant states and T(j — > l + l~) 
denote the partial decay width for the transition of vector charmonium state to massless 
lepton pair, which can be expressed in terms of the decay constant of charmonium through 
the relation [76] 



r(j -> £ + r) 



16 /. 



2 



ira 



2 -iu Jj 

em 27 Mj 



The phenomenological parameter kj in Eq.(2.5) is to account for inadequacies of the fac- 
torization approximation, and it can be determined from 

BK{B C D* s J/ij) D* s £ + £~) = BK{B C £>* J/V>) ■ BTZ(J/ip £ + £~). 

The function iOj(q 2 ) introduced in Eq.(2.5) is to compensate the naive treatment of long 
distance contributions due to the charm quark loop in the spirit of quark-hadron duality, 
which can overestimate the genuine effect of the charm quark at small q 2 remarkably . 
The quantity ujj(q 2 ) can be normalized to uij(M^) = 1, but its exact form is unknown 
at present. Since the dominant contribution of the resonances is in the vicinity of the 
intermediate ipi masses, we will simply use ujj(q 2 ) = 1 in our numerical calculations. 

Moreover, the non factorizable effects from the charm quark loop brings further cor- 
rections to the radiative transition b — > 57, and these can be absorbed into the effective 
Wilson coefficients Cy which then takes the form [76, 82-861 



c 7 e// 0u) = c 7 (/i) + cu S7 Gu) 



with 



lv U/23 (Gi(x t ) - 0.1687) - 0.03C 2 (aO 
y 



(2.7) 



x t (x 2 - 5x t - 2) 3x 2 In 2 x t 
<{x t -l? + A(x t -1) 4 



Gl (x t ) = ^ ^1 + ^i_iL (2 . 



where r\ = a s {my/)la. s {\i), xt = m 2 /m 2 v and Cb^-sy is the absorptive part for the b — > 
sec —> 57 rescattering. 

2.2 Weak Annihilation Amplitude 

The charmed B-meson (B c ) is made up of two different heavy flavors, 6-quark and c-quark, 
which brings WA contributions into the play. Using the procedure developed in refs. [87, 88] 
for B c — > .D*7, the WA amplitude for the decay B c — > D*£ + £~ can be written as 

M^ DH+e - = ^V cb V c * s \-ie^e*y a qP)Cr n (q 2 ) + (e ■ qp, + p ■ qe,) IC a 2 nn (q 2 )} l-fl 



2y/2- 



7T 



(2.9) 

where /C" nn (g 2 ) and K,2 nn (q 2 ) are the weak annihilation form factors. 



1 For a more detailed discussion on long-distance and short-distance contributions from the charm loop, 
one can refer to references [9, 76-82]. 
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Before proceeding further we would like to mention that we use parametrizations for 

i.e. 



weak annihilation form factors /C" nn (g 2 ) and K.f in {q 2 ). 

/q™(0) 



yanni l 



l + a 



(2.10) 



The values of JCf nn (0) and K,?, nn {0) are calculated by using QCD sum rules [37, 38] and the 
values of the parameter a and /3 are given in Ref. [26-31, 37, 38], which are summarized 
in Table 1. 



K ann (q 2 ) 


V~ann 
^0 


a 




fcann ( g 2j 


0.23 


-1.25 


-0.097 


K,% nn {q 2 ) 


0.25 


-0.10 


-0.097 



Table 1. B c — >• D* form factors corresponding to WA in the QCD Sum Rules. /CJ™(0) denote the 
value of form factors at q 2 = while a and /3 are the parameters in the parametrizations shown in 
Eq. (2.10) [38]. 



2.3 Parameterizations of the Matrix Elements and Form Factors 

The exclusive B c — > D*£ + £~ decay involves the hadronic matrix elements which can be 
obtained by sandwiching the quark level operators give in Eq. (2.3) between initial state 
B c meson and final state D* meson. These can be parametrized in terms of form factors 
which are the scalar functions of the square of the four momentum transfer (q 2 = (p — k) 2 ). 
The non vanishing matrix elements for the process B c — > D* can be parametrized in terms 
of the seven form factors as follows 

(D* s (k,e) \s lfl b\ B c (p)) = Mz V +M D * W^V^ (2.11) 

(D* s (k, e) \s W b\ B c (p)) = i (M B - + M DJ ) e*»A (q 2 ) - i J ^0_ ^* . p){p + k f 

A_ ( a 2) 

(e*-p)q» (2.12) 



M Bc + M D * 



(D*(k,s)\s(l± l5 )b\B c ( P )) = T i 2Md ° (e*-p)A (q 2 ) (2.13) 

m b + m s 

where p is the momentum of the B c , e and k are the polarization vector and momentum 
of the final state D* vector meson. Whereas A${q 2 ) can be parametrized as 

In addition to the above form factors there are some penguin form factors, which we 
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can write as 



(D* s (k,e)\sa tlu q"b\B c (p)) = 2iT 1 (q 2 )e flua f ) e* v p a kf ! 



(D* s (k,e)\sa fiU q" 1 5 b\B c (p)) 



M 



Be 



M 



D* 



(2.15) 
(2.16) 



+(e*-p) 



Ml 



Uq 2 ). 



The form factors Ay (q 2 ) , A (q 2 ) , A + (q 2 ) , A. (q 2 ) , T x (q 2 ) , T 2 (q 2 ) , T 3 (q 2 ) are the 
non-perturbative quantities and to calculate them one has to rely on some non-perturbative 
approaches and in our numerical analysis we use the form factors calculated by using QCD 
Sum Rules [38]. The dependence of these form factors on square of the momentum transfer 
(q 2 ) can be written as 



l + a- 



(2.17) 



'Or 



"To 



where the values of the parameters F (0), a and b is given in Table 2. 



F(q 2 ) 


F(0) 


a 


b 


Ay (q 2 ) 


0.54 ±0.018 


-1.28 


-0.23 


A (q 2 ) 


0.30 ±0.017 


-0.13 


-0.18 


A + (q 2 ) 


0.36 ±0.013 


-0.67 


-0.066 


A-(Q 2 ) 


-0.57 ±0.04 


-1.11 


-0.14 


Ti(q 2 ) 


0.31 ±0.017 


-1.28 


-0.23 


T 2 (q 2 ) 


0.33 ±0.016 


-0.10 


-0.097 


T 3 (q 2 ) 


0.29 ±0.034 


-0.91 


0.007 



Table 2. B c — > D* form factors corresponding to penguin contributions in the QCD Sum Rules. 
F(0) denotes the value of form factors at q 2 = while a and b are the parameters in the paramctriza- 
tions shown in Eq. (2.17) [38]. 

Prom Eq. (2.3) it is straightforward to write the penguin amplitude 
<peng _ Gfu 



v tb v t * s [Ti(h»e) + T 2 {h^H) + r 



(2.18) 



where 



Tl = fi(q 2 )e^ a pe*»p a kP-if 2 (q 2 )e; + if 3 (q 2 )(e* -p)P^ (2.19) 
T 2 = h{q 2 )^e™p a k? -ih{q 2 )el + if 6 ( q 2 )( £ * ■ p)P )1 + if 7 ( q 2 )(e* ■ p)P fl (2.20) 
T =2ih(q 2 ){e* -p) (2.21) 

with P^=p^ + k^. 

The next task is to calculate the decay rate and the helicity fractions of D* meson 
in terms of these auxiliary functions which contains both long distance (form factors) and 
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short distance (Wilson coefficients) effects and these can be written as 



/i(. 2 ) =4(C?" + C^f-^^T^) + (C^ + C'^) ^l q 2 Dt (2-22) 



/ 2 (g 2 ) = 2 (c?" - c?")^^i T ^) (m! c - M 2 D , 

+ (C e 9 ff - C' 9 eff )A (q 2 ) (M Bc + M D *) (2.23) 

f 3 ( q 2 ) =4(C°/f - C*")=*^ Ltf) + q 2 , T3(g2) 

9 V ( M 1 c - M d,) 



M Bc + Md» 



+ (C 9 e// - C' 9 etf ) „ + , (2.24) 



/ 5 (9 2 ) =2(Cio - Cj„)4)(? 2 ) {M Bc + M D .) (2.26) 
Atf) =2( Cl0 - C; 0>S ^; (2.27) 

M 9 2 ) =4(C 10 - CW^^; + - cy-J^-j^tf) (2.28) 

*(«•) = - (09. -C&J^^jAtf) (2.29) 
Here the NHBs contribution are encoded in the auxiliary functions fa and f$. 

3 Physical Observables for B c -> D*£ + £~ 

In this section we will present the calculations of the physical observables such as the 
branching ratios BIZ, the forward-backward asymmetries Afb, the lepton polarization 
asymmetries Pl,n.t and the helicity fractions /tj of D* meson in B c — > D*l + l~ decay by 
incorporating both the weak annihilation (WA) and the penguin amplitudes. 

3.1 The Differential Decay Rate of B c -> D* s t + £~ 

In the rest frame of B c meson the differential decay width of B c — > D*£ + £~ can be written 



as 



dT{B c ^ D*J+£~) 1 1 



dq 2 " (2vr) 3 32M| c J_ u{q2) 

where 



du\M\ 2 (3.1) 



M = M WA + .M PENG (3.2) 
q 2 = (PI++PI-? (3.3) 
u = (p-Pi-f - (p-Pi+f (3.4) 
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Now the limits on q 2 and u are 



4m 2 < q 2 < (Mb c — Mfl») 2 
-u(q 2 ) <u< u(q 2 ) 



with 



and 



u{q 2 



A 1 



4m 2 \ 



(3.5) 
(3.6) 

(3.7) 



A = \{Ml c ,Ml t , q 2 ) = M% c + M% t + g 4 - 2M 2 B M 2 Dt - 2M 2 D ,q 2 - 2q 2 M 2 Bc 

Here m corresponds to the mass of the lepton which for our case are the \i and r. The 
total decay rate for the decay B c — > D*£ + £~ can be expressed in terms of WA, penguin 
amplitude and the interference of these two, which takes the form 



dT _ dr^A q(p PENG dr WA-PENG 
dq 2 dq 2 dq 2 dq 2 



with 



dr WA 

dq 2 
^p PENG 



G 2 F \V tb V t */a 2 



2 u TT 5 3M Br M 2 ) ,q 2 



u(q 2 )xA{q 2 ) 



G 2 F \VtbVts\ 2 a 2 ( 2) B , 2} 



dT WA- 



PENG 



G 2 F \V cb V* s \\V tb V t * s \a 2 
2 u n 5 3M Bc M 2 )t q 2 



dq 2 



■u(q 2 ) x C (q 2 ) 



(3.8) 

(3.9) 
(3.10) 
(3.11) 



The function u(q 2 ) is defined Eq. (3.7) and A(q 2 ), B{q 2 ) and C(q 2 ) are 



A (q 2 ) = i (2m 2 + q 2 ) k 2 



%\M 2 D *q 2 (/C? nn ) 2 + m nn ) z { V2M 2 D *q l (\ + \M 2 B q 



2 „2/ 



r2 „2^ 



+ A z + A ( A + 4q 2 M 2 D * + 4g 



(3.12) 



B(q') = 8M 2 D ,q 2 Xi (2m' + q')\h(q 



,2\l 2 



(4m 2 - g 2 ) 



1/4(9 



2M 2 



+ 4Mf)*q 2 i [2m 2 +q 2 ) 



x (3 |/ 2 (g 2 )| 2 - A |/ 3 (g 2 )| 2 ) - (4m 2 - q 2 ) ( 3 |/ 5 ( 9 2 )| 2 - A |/ 6 (g 2 )| 2 ) 



+ \(2m 2 +q 2 ) f2(q 2 ) + [M 2 Bc -M 2 D ,-q 2 )f 3 (q 2 ) + 24m 2 M 2 D ,X\f 7 (q 2 ) 



-2M 2 



- (4m 2 - q 2 ) h{q 2 ) + (M 2 c - M 2 . - q 2 ) f 6 (q 2 ) + (g 2 - 4m 2 ) A |/ 8 (g 



JM2 



12mV [K(/b/ 7 *)-K(/6/ 7 *)] 



C(g 2 ) = 2 K / 2 (g 2 )/Cr <? 2 (2m 2 + (? 2 )|a + 6M 2 , (m 2 c - M 2 . + g 2 ) | 

- kA [2/ 1 (g 2 )/Cr n M 2 »g 4 + f 3 (q 2 )IC a 2 nn (2m 2 + g 2 ) (A + q 4 + 4M Bc M D *) 



(3.13) 



. (3.14) 



- 11 - 



where 

_ ^ 2 M Dt f B J Dt 

K — — 5 o\ 2 ' (^O.IOJ 

(nic — iTig)q z 
3.2 Forward-Backward Asymmetries 

The differential forward-backward asymmetry Afb of final state lepton for the said decay 
can be written as 

dA FB {s) f 1 d 2 r j a r° d 2 r 



dq 2 



r 1 d z r r d 2 r 

Jq dq 2 dcos9 J_idq 2 dcos9 ^ 



From experimental point of view the normalized forward-backward asymmetry is more 
useful, defined as 

r 1 f/ j cose - f\ . £ T . dcosfl 

. JO dq z d cos J — 1 dq z d cos 

ApB = fl 3 a 

i j -L a d cos 8 

J — 1 ogr a cos o 

The normalized Afb fo r -Sc — * D*£ + £~ can be obtained from Eq. (3.1), as 
•Afb = ^ L o s I WJ 2 g\(g 2 )(2^(/ 4 /Cr rt )(Ml e - M 2 + g 2 ) + 2/cM(/ 5 /Cr n ) 



dr/dg 2 2 11 tt 5 to| 

+4i?e[/ 2 7 4 + /r/s] + 2AK[/ 3 7 8 ] + 4K[/ 2 /|] (-M 2 c + M 2 , + g 2 ' 
+2 K K[/ 8 ]/Cr n (A + (M 2 + q 2 ) M 2 c - q 4 ) ) (3.17) 



where k is defined in Eq. (3.15) and dT/dq 2 is given in Eq. (3.8). 
3.3 Lepton Polarization Asymmetries 

In the rest frame of the lepton i~ , the unit vectors along longitudinal, normal and transver- 
sal component of the t~ can be defined as [47, 48, 89]: 

s -^ = (0,^)=(0,J^Y (3.18a) 
s -» = (0, gp) = I 0, A^l— I ; (3.18b) 



T 

k x p. 

s/ = (0, e- N ) = (0, e r x e L ) , (3.18c) 

where p_ and are the three-momenta of the lepton £~ and D* meson respectively in the 
center mass (cm.) frame of system. Lorentz transformation is used to boost the 

longitudinal component of the lepton polarization to the cm. frame of the lepton pair as 

/ CM \ m m \p_ | / 

where E and m are the energy and mass of the lepton. The normal and transverse com- 
ponents remain unchanged under the Lorentz boost. The longitudinal (Pl), normal (Pn) 
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and transverse (Pt) polarizations of lepton can be denned as: 

dr 

dq- 



dq 



(3.20) 



where i = L, N, T and is the spin direction along the leptons P. The differential 
decay rate for polarized lepton P in B c — > D*£ + i~ decay along any spin direction ^ is 
related to the unpolarized decay rate (3.8) with the following relation 



dq 2 



1 ( dT 



2 \dq< 



1 + + P%(% + pTe$) • p 



(3.21) 



The expressions of the numerator of longitudinal, normal and transverse lepton polariza- 
tions can be written as 



P L (q 2 ) oc 



4A 



4m 2 



3M 2 



x u(Ml) i + 



12g 2 M 2 

^ s 

A 



+ A3?(/ 3 /e) + 8g 2 Af|, ; K(/i/| 



-M 2 c + M^, + g 2 J [5R(/ 3 /|) + R(/ 2 /*)] 

K(/ 5 /|) + »(/ 6 /|) (-M| c + M|,j) - q 2 K(hfz) 
+ AMl*KR(1Ci nn fl) + [A + g 2 (AfJ, + M 2 D , - g 2 )] KM{K% nn ft 
6g 2 M 2 



3 

H — m 
2 



+ 



M 



M 



A 



K K(/cr n /2) 



(3.22) 



-M 2 +M 2 +g 2 



q 2 W2fj) + (M 2 c - M 2 D . - g 2 )K(/ 2 / 5 *) + (g2 2 ^ m2) ^(/ 5 /|; 



+ Zq 2 M 2 D M{hr 2 ) + A( g 2 _ 4m 2 )5ft(/ 6 / 8 *) + 



A 



- 2 / 3M|>, + M| + - 2 



4g 2 M 2 j^(X:r n /|) - 4g 2 M 2 ; (M 2 c - M 2 * + g^^/i*; 
A + g 2 (M 2 c + ML - r 2 ' 



« ( W"7 5 *) - (M 2 c - M 2 ,)K(/Cr n /e) + g 2 K(/Cr7 7 *) 
+ 4K 2 g 2 Af|,, (m% c - M 2 D , + g 2 ) /Cr n /Cr n | (3.23) 



miT\/ (q 2 — 4m?) A 

^ (g 2 ) oci L_ ^ M 2 , [43(/ 2 / 4 *) + 49(/ 1 /*) + 2^(/ 5 /Cr n )] 



M 2 



, N (A + g 2 (M 2 +M 2 » -g 2 )) 
+ 2M 2 «3f(/ 4 /Cr n ) (Mi c - M 2 + g 2 ) - A3(/ 6 / 7 *) - ^ 1- U^{KT n ft) 



2m 



+ (-M 2 c + M 2 D . + g 2 



»(/7/l) + a^^Ws) 



+ (M 1 c + 3M 2 . - g 2 ) 3?(/ 5 / 6 *) + ^{hft) 



(3.24) 
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where k is defined in Eq. (3.15) along with auxiliary functions /i, /2, • " " ; fs m Eqs. (2.22- 
2.29). Here we have dropped out the constant factors which are however understood. 



3.4 Helicity Fractions of D* in B c D*£ + £~ 

We now discuss helicity fractions of D* in B c — > D*£ + £~ which are interesting variable and 
are as such independent of the uncertainties arising due to form factors and other input 
parameters. The final state meson helicity fractions were already discussed in literature for 
B —7- K* (Ki) £ + £~ decays [90-93]. Even for the K* vector meson, the longitudinal helicity 
fraction fi has been measured by Babar collaboration for the decay B — > K*£ + £~ (Z = e, /u) 
in two bins of momentum transfer and the results are [94] 



h 



0.77toH ± °- 07 ' 



f L = 0.51±g;l| ± 0.08, 



0.1 < q 2 < 8AlGeV 2 
q 2 > W.2AGeV 2 



(3.25) 



while the average value of fi in full q 2 range is 



f L = 0.631^ ± 0.05, q z > 0.1GeV~ 

The explicit expression of the decay rate for B~ — > D*~£ + - 
terms of longitudinal Tl and transverse components 1~Y as 



(3.26) 

decay can be written in 



dF L (q 2 ) _ dT™ A (q 2 ) dr™ G (g 2 ) ^pWA 



dq 2 
dT ± (q 2 ) 
dq 2 

dT T (q 2 ) 
dq 2 



PENG/^2 



dq 2 



+ 



dq 2 



drf A ( q 2 ) drl™ G (g 2 ) 

dq 2 
dT + (q 2 ) 



dq 2 



+ 



dq 2 
dT.(q 2 ) 
dq 2 ' 



+ 



+ 



dq 2 
^WA-PENG 



dq 2 



(3.27) 
(3.28) 
(3.29) 



where 



dq 2 
d pPENG ((? 2 ) 



dq 2 



dT WA- 



PENG/^2 



(q 2 ) 



dq 2 

dr¥ A (q 2 ) 

dq 2 
dTl™ G (q 2 ) 



dq 2 

^pWA-PENG 



i<f) 



dq 2 



G 2 F 



G 2 f 



G\ 



G 2 f 



V cb V* s \ 2 a 2 u{q 2 ) 1 



2 11 vr 5 



V tb V t * s \ 2 a 2 u{q 



2 U 7T 5 
cb^ cs 



M l 



x- 3 A L 



WA 



X ~ A L 



PENG 



V rh V*\ \V tb V t *\a 2 u{q 2 ) 



2 11 vr 5 
V cb V c *fa 2 u(q 



Ml 



X I^WA-PENG 



2 11 7T 5 



Ml 



V tb V t * s \ 2 a 2 u(q 2 



2 
3" 

4 



Ml 



x- 3 A ± 



PENG 



2 n 7r 5 

V cb V* s \\V tb V t * s \a 2 u(q 2 



2 n 7r 5 



Ml 



2 
3" 



X ^A-EP. 



(3.30) 
(3.31) 
(3.32) 
(3.33) 
(3.34) 
(3.35) 
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The different functions appearing in Eqs. (3.30-3.35) can be expressed in terms of 
auxiliary functions (cf. Eqs. (2.22-2.29)) as 



I (£i nn (q 2 )) |? 2 A(A + Aq 2 Ml t ) - 4M 2 A(2A + 8q 2 M 2 D *) 

D s 

- q 2 (Ml - M 2 , - q 2 )' (A - 2uV)) } + (JCT tf)f { ~ AV ~ ^) 
+ 12Xq 2 ((M 2 Bc - M 2 D ,) 2 - M 2 Dt ) + q 2 (8q 2 M 2 D , - \){M 2 Bc - M 2 D . + q 2 ) 2 

- 2u 2 (q 2 ) q 2 ((M 2 Bc - M 2 ) 2 + q') + 4m 2 ((M 2 c - M 2 ) 2 - Ml (3.36) 



^PENG 



q 2 M 2 



24 \f 7 (q 2 )\ 2 m 2 M 2 D . A + (2m 2 + q 2 ) (M 2 Bc - M 2 - q 2 )f 2 (q 2 ) + A/ 3 (g 2 ; 



+ (q 2 - 4m 2 ) (M 2 - M 2 - q 2 )h{q 2 ) + A/ 6 ( <? 2 



+ i(g 2 -4m 2 )A|/ 8 | 2 



(3.37) 



^WA-PENG 



q 2 M 2 



Re(h (q 2 ) KT n {q 2 ))\ (A + 4M 2 ,g 2 ) (8m 2 V\ + q 2 (2u(q 2 ) 



- 4M 2 ^ 2 A j + Re(f 2 {q 2 ) KT n {q 2 ))U 2 u 2 {q 2 )(M 2 Bc - M 2 D , - q 2 ) 

+ Qq 2 \{M 2 D , - M 2 Bc ) + q 2 (X - 8g 2 M 2 *)(M 2 c - M 2 , + q 2 ) - Am 2 q 2 (4q 2 M 2 D , + A)| 

+ Re(f 3 {q 2 ) /CrV)){A 2 (4m 2 - q 2 ) + q\q 2 u{q 2 )^f\ - 6A(M 2 c + M 2 ,)) 



+ q\M B -Mf n )^X-u 2 {q 2 )) 



iWA 



^PENG 



(2m 2 + q 2 ) 



A(/Cf"V)) 2 + {lC% nn (q 2 )y (A + 4M£.g 
2 



,2N\ 2 



2 „2 



(3.38) 
(3.39) 



i ± -(g 2 -4m 2 ) f 5 (q 2 )TVXh(q 2 ) + {q 2 + 2m 2 ) f 2 (q 2 ) ± VA/i(<? 2 ) (3.40) 
A WA -PENG = _ J 2 ,/X (g 2 - 4m 2 ) J Re(/ 2 (g 2 ) £rV)) + 4A(g 2 + 2m 2 )Re{h {q 2 ) KfV)) 

± 2(</ 2 + 2m 2 )(M 2 c - M 2 , + g 2 )[2i?e[(/l (g 2 ) JCfV))] ^A T 2i2e(/ 2 (g 2 ) /Cf»V))] 

(3.41) 

Finally the longitudinal and transverse helicity amplitude becomes 



f±(q 2 ) 



_ dT L {q 2 )/dq 2 



dT(q 2 )/dq 2 
dV ± {q 2 )/dq 2 



dT{q 2 )/dq 2 
h(q 2 ) = f+(q 2 ) + f-(q 2 ) 



(3.42) 
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so that the sum of the longitudinal and transverse helicity amplitudes is equal to one i.e. 
/l(^ 2 ) + friq 2 ) = 1 for each value of g 2 [90]. 

As we have disscussed in the Introduction that the WA contribution plays a cruicial 
role in the B c — > D*£ + £~. This feature can also be seen via the expressions of different 
observables for example branching ratio which is given in eq. (29), there are three terms 
A(q 2 ), B(q 2 ) and C(q 2 ) correspond to WA, Penguin and the cross term of Penguin and 
WA, although the WA contribution A(q 2 ) trying to suppress the effects of NP but the 
term C(q 2 ) some how compensate and enhanced the NP effects. Same is the case for the 
other calculated observables. Therefore, the interference term between the NP and WA 
play a cruicial role and make the NP effects distinct form the SM. 

4 Numerical Results and Discussion 

We present here our numerical results of the branching ratio (BIZ) , the forward-backward 
asymmetry (Afb), lepton polarizations asymmetries Pl,n,t and the helicity fractions 
(/l,t) of D* for the B c — > D*£ + £~ decays with £ = fi,r. The numerical values of the 
input parameters which are used in the subsequent analysis are summarized in Table 3: 

m Bc = 6.277 GeV, m D . = 2.1123 GeV, m b = 4.28 GeV, = 0.105 GeV 

m T = 1.77 GeV, \V cb V* s \ = 4.15 x 10" 2 , \V tb V t * s \ = 4.1 x 10~ 2 , t b = 0.453 x 10" 12 sec 

G F = 1.17 x 10~ 5 GeV" 2 , a" 1 = 137, f B = 0.35 GeV, f D = 0.30 GeV 

Table 3. Values of input parameters used in our numerical analysis [95]. 



Wilson coefficients 




r eff 


Cio 


SM 


-0.313 


4.334 


-4.669 


SUSY I 


+0.376 


4.767 


-3.735 


SUSY II 


+0.376 


4.767 


-3.735 


SUSY III 


-0.376 


4.767 


-3.735 


SUSY SO(10) (A = -1000) 


-0.219 


4.275 


-4.732 



Table 4. Wilson Coefficients in SM and different SUSY models but without NHBs contributions 
[45-48]. 



Wilson coefficients 






°10 


SM or SUSY I,II,III 











SUSY SO(10) (A = -1000) 


0.039 + 0.038i 


0.011 + 0.072i 


-0.075 - 0.67i 



Table 5. Primed Wilson Coefficients in SM and different SUSY models but without NHBs contri- 
butions. Where the primed Wilson coefficients corresponds to the operators which are opposite in 
helicities from those of the SM operators. 

The values of Wilson coefficients used in our numerical analysis are taken from Refs. 
[45-48] and summarized in Tables 4, 5 and 6. In the following analysis, we will focus on 
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Wilson coeff. 






c q 2 




SM 














SUSY I 














SUSY II 


6.5(16.5) 





-6.5(-16.5) 





SUSY III 


1.2(4.5) 





-1.2(-4.5) 





SUSY SO(10) 


0.106 + 0i 


-0.247 + 0.242? 


-0.107 + 0i 


-0.25 + 0.246i 


(A Q = -1000) 


(1.775 + 0.002i) 


(-4.148 + 4.074i) 


(-1.797 -0.002?) 


(-4.202 + 4.128i) 



Table 6. Wilson coefficient corresponding to NHBs contributions in different SUSY scenarios. 
Where the primed Wilson coefficients are for the primed operators from NHBs contribution in 
SUSY SO(10) GUT model. The values in the parentheses are for the r lepton case. 



the parameter space of large tan /3, where the NHBs effects are important owing to the fact 
that the Wilson coefficients corresponding to NHBs are proportional to (m^m/m^) tan 3 /3, 
with h = h°,A°. Here, the tan/3 contributes from the chargino-up-type squark loop and 
the tan 2 /3 appears from the exchange of the NHBs. In the Ref. [62, 63] it is pointed out 
that at large value of tan (3 the Cq, compete with and can outpace C^' in some region. 
Depending on the magnitude and sign of the SUSY parameters one can think of many 
options in the parameter space, but experimental results i.e., the decay rate of b — > S7 and 
b — > s£ + £~ restrict us to consider the following scenarios for MSSM: 

• SUSY I: refers to the regions where SUSY destructively contributes and changes 
the sign of C% , which will have drastic effects on the observables, but without 
contribution of NHBs. 

• SUSY II: corresponds to the region where tan f3 is large and the masses of the spar- 
ticles are relative small. 

• SUSY III: points to the regions where tan /3 is large and the masses of superpartners 
are also relatively large, i.e. > 450 GeV 2 . 

In SUSY I and SUSY II scenarios one can accommodate the non-zero crossing of the 
forward-backward asymmetry in B — > K* ji + fi~ decay [52-55]. Since the primed Wilson 
coefficients are for the primed operators in Eq. (2.2) which appear in the SUSY SO(10) 
GUT model. In Table 6 it is mentioned that the values of Cq_ is different for the different 
choice of final state lepton, which is due to the fact that contributions from the NHBs are 
proportional to the lepton mass. 

The other two conditions, apart from the large tan f3 limit, which are responsible for 
the dominant contributions from NHBs are: (i) the mass values of the lightest chargino 
and lightest stop should not be too large and (ii) the mass splitting of charginos and 
stops should be relatively large, which also indicate large mixing between stop sector and 
chargino sector as pointed out in Ref. [45]. Since the SUSY effects are very sensitive to 
the sign of the Higgs mass term and SUSY contributes destructively appears when the 
sign of this term becomes minus. In Ref. [45] it is shown that there exist some regions 
of SUSY parameter space where NHBs can dominantly contribute to the FCNC process 
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b — > s£ + £~ due to flip of the sign of C7 from positive to negative, within the constraint on 
b 57. Also it is noticed that when the masses of sparticles are relatively large, say about 
450 GeV, there exits considerable regions in the SUSY parameter space where NHBs could 
contribute dominantly as in the case of SUSY III and the SUSY SO(10) models. However, 
in these scenarios the sign of Cj remains unaltered with respect to the SM sign because of 
the cancellation of the contributions of charged Higgs and charginos with each other. 

It is worth mentioning here that the decay B q — > £ + £~ is a clean channel to probe for 
the NHBs effects in SUSY models at large tan/3. Since its branching ratio in the SUSY 
models can be written as [96] 



where Cq 1 and Cq 2 are referring to the NHBs contributions which are absent in the SM but 
in MSSM they are proportional to {m^m/rn^) tan 3 /?, with h = h°, A . The value of C\q is 
large in the SM but it is suppressed by the factor 2m/mB s , therefore, the corresponding 
SM value of branching ratio is [97] 



The branching ratio corresponding to different values of the NHBs parameters can be 
enhanced by a factor ranging from 10 1 — 10 3 , but the recent upper bound by CDF col- 
laborations, 3.9 x 10 -8 at 95% confidence level [97], is about 12 times larger than that of 
SM prediction. It can be noticed from Eq. (4.1) that the branching ratio for B s — > 
is directly proportional to the NHBs contributions. Therefore, this stringently constraints 
the parameter space of the SUSY models and especially the large value of the Wilson Coef- 
ficients Cq 1 and Cq 2 are severely constrained. Therefore, the precise measurement of this 
decay at the future experiments will help us to get useful constraints on the SUSY param- 
eters and if considerably large deviation from SM predication is measured, then along with 
the signal of supersymmetry, this would have important implications on the Higgs searches 
at LHC. 

The purpose of the present study of the SUSY effects in B c — > D*£ + £~ , with £ = fi,r, 
is to incorporate the constraints provided by B s — > as well as to see the effects of 

NHBs at the larger extent in these FCNC processes. Since this decay mode is unique in 
a sense that it contains WA contributions along with the penguin and so it is interesting 
to see how SUSY affects different physical observables in this process. For our numerical 
analysis we have set the values of Cq 1 and Cq 2 between to 6.5(16.5) for muons and 
(tauons) to check the dependency of different observables on these NHBs contributions. 
We hope that in the future experiments such as the Tevaran and the LHC, with the more 
data on these decays, will help us to test more precisely the constraints obtained from the 
decay channel B s — > fi~ . 




(4.1) 




(4.2) 
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First, we discuss the branching ratios (BIZ) for the decays B c — > D*£ + £~ , with £ = fi,r, 
which we have plotted as a function of q 2 (GeV 2 ) in Fig. 1, without and with long-distance 
contributions in the Wilson coefficients, both in the SM and in the SUSY scenarios. This 
figure depict that the values of BIZs, both for the case of muons and tauons as final state 
leptons, get sizeably influenced due to the SUSY effects which come through the new 
parameters, i.e the modified and new Wilson coefficients corresponding to the operators 
described in Eq. (2.2). One can see clearly from these graphs that the increment in the 
values of BTZs in the SUSY I model is due to the relative change of the sign of Cf ff with 
respect to that of the SM; while the large deviation in SUSY II model from the SM values 
is mainly due to the contributions of NHBs and due to the relative change of the sign of 

. Since the values of Wilson Coefficients corresponding to NHBs is small for the cases 
of SUSY III and SUSY SO(10) GUT models, so one should expect small deviation in the 
BTZs from the SM values and this signature is clear in Fig. 1. Moreover, the NP effects due 
to the SUSY, manifest in the BTZs throughout the q 2 region irrespective of the mass of the 
final state leptons. In addition, one can also extract the constructive behavior of SUSY I 
and SUSY II to the BTZ from Table 7. Furthermore, we have also plotted the BTZs with 
the long-distance contribution corresponding to the SM and SUSY models in Fig. l(b,d). 



(a) (b) 




13 14 15 16 17 13 14 15 16 17 

? 2 (GeV 2 ) 9 3 (GiiV 2 ) 



Figure 1. The dependence of branching ratio of B c — > D* s £ + £ (£ = (m, t) on q 2 without long- 
distance contributions (a,c) and with long-distance contributions (b,d) for different scenarios of 
MSSM and SUSY SO(10) GUT model. In all the graphs, the solid, dashed, dashed-dot, dashed- 
double dot and dashed-triple dot curves correspond to the SM, SUSY I, SUSY II, SUSY III and 
SUSY SO(10) GUT model, respectively. 
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Branching 


B c D* s ^ 


> 


B c -> 






"r 


B c -> £>>+ 


/' 


ratios 


Without LD 


With LD 


Without LD 


With LD 


SM 


1.60 x 10' 


-6 


2.81 x 10' 


-5 


1.33 x 10- 


-8 


1.19 x 10- 


-6 


SUSY I 


1.91 x 10' 


-6 


2.85 x 10" 


-5 


3.06 x 10- 


-8 


1.21 x 10- 


-6 


SUSY II 


2.29 x 10" 


-6 


2.89 x 10" 


-5 


2.89 x 10- 


-8 


1.22 x 10- 


-6 


SUSY III 


1.58 x 10" 


-6 


2.81 x 10" 


-5 


1.35 x 10- 


-8 


1.19 x 10- 


-6 


SUSY SO(10) 


1.60 x 10" 


-6 


2.81 x 10- 


-5 


1.52 x 10- 


-8 


1.19 x 10- 


-6 



Table 7. Branching ratio for B c -> D*£ + £~ (£ = ( u, r) in the SM and different SUSY models. 



It is important to mention that as an exclusive decay, there are different sources of 
uncertainties involved in the analysis of the above mentioned decay. The major source 
of uncertainties in the numerical analysis of B c — > D*£ + £~ {£ = fi, r) decays originate 
from the B c — > D* transition form factors calculated in the QCD sum rule approach 
[38] as summarized in Table 2. But it is also important to stress that these hadronic 
uncertainties have almost no influence on the various asymmetries including the forward- 
backward asymmetries, lepton polarization asymmetries and helicity fractions of D* in the 
decays B c — > D*£ + £~ because of the cancellation among different polarization states and 
this make them a good tool to probe beyond the SM. 

To illustrate the generic effects due to the SUSY models on the forward-backward 
asymmetry Afb, we plot as a function of q 2 in Fig. 2. For the zero position of Afb 

it is also argued that the uncertainty in the zero position of the Afb due to the hadronic 
uncertainties is negligible [9] to leading order in a s . Therefore, the zero position of the 
Afb can serve as a stringent test for the NP effects including the SUSY. Figures 2(a) and 
2(b) describe the Afb for B c —> D*[i + n~~ with and without long-distance contributions in 
the Wilson coefficients respectively, where the different SUSY models show clear deviation 
from the SM curve. Figure 2(a) clearly shows that SUSY I and SUSY II do not cross the 
zero of the Afb because of the fact that the zero crossing of Afb of SM is solely due to the 
opposite signs of C% ff and Cg ff but in SUSY I and SUSY II the signs of these coefficients 
are same and hence the zero point of Afb disappears in both cases. Furthermore, in the 
SUSY III model, due to the opposite sign of C% and Cg^, the Afb passes from the 
zero but the zero position shifts mildly to the right (about 0.2 GeV 2 ) from that of the SM 
value, i.e 3.3 GeV 2 , due to the contribution from the NHBs. Moreover, the SUSY SO(10) 
GUT model also shows the similar behavior except that the zero position of the Afb shifts 
considerably to the left (about 1 GeV 2 ) from that of the SM value. Hence, the precise 
measurement of the zero position of Afb f° r the decay B c — > will be a very good 

observable to yield any indirect imprints of NP due to SUSY models and can serve as a 
good tool to distinguish among the different variants of SUSY models. 

Besides the zero position of Afb, the magnitude of Afb is also an important tool to 
establish NP. Particularly, this con be more exciting when the tauons are the final state 
leptons, where the zero of the Afb is absent. In the Figs. 2(c) and 2(d) the Afb of 
B c —7- D*t + t~ is plotted with and without long-distance contributions, respectively, where 
the different SUSY models present a distinct variations from that of the SM magnitude. 
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Figure 2. The dependence of forward-backward asymmetry for the decays B c — > D*£ + £ (£ = 
/i, r) on q 2 without long-distance contributions (a,c) and with long-distance contributions (b,d) for 
different scenarios of MSSM and SUSY SO(10) GUT model. 



A closer look on the pattern of Fig. 2(c) indicates that the SUSY I and SUSY II models 
decrease the magnitude of Afb from its SM value. Whereas, SUSY III has very mild 
deviation from that of the SM behavior but SUSY SO (10) GUT model has an increment in 
the magnitude of Afb compared to that of the SM value for the case where tauons are in 
the final state. It is valuable to comment here that just like the zero position of the Afb, 
the magnitude of Afb depends on the values of the Wilson coefficients , Cg , Cio and 
Cq. and the effects due to the hadronic uncertainties on the magnitude of Afb are almost 
insensitive. Hence, the deviation in the magnitude of Afb due to the SUSY parameters 
is prominent and can be measured at the experiments which indeed help us to understand 
the constraints on the parameter space of the SUSY. 

We now consider another interesting observable to get the complementary informa- 
tion about NP in B c — > D*£ + £~ {£ = /i, r) decays, i.e. the lepton polarization asymmetries 
which are shown in Figs 3, 4 and 5. Since different polarization asymmetries depend on the 
Wilson Coefficients so one can expect large dependency of these asymmetries on different 
SUSY variants and hence making these observables fertile to extract the NP. Therefore, 
we expect that even in the SUSY I model where the value of the Wilson coefficient corre- 
sponding to NHBs is zero, the values of these polarization asymmetries would be mainly 
modified from those of the SM value, because of the change in the sign of the terms pro- 
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Figure 3. The dependence the probabilities of the longitudinal lepton polarization asymmetries, 
Pl, for the decays B c — > D*£ + £~ [£ = /i, t) on q 2 without long-distance contributions (a,c) and with 
long-distance contributions (b,d), respectively, for different scenarios of MSSM and SUSY SO(10) 
GUT model. 



portional to Cj C\q. Now we focus on the longitudinal polarization asymmetry for the 
decays B c — > D*£ + £~ (£ = fi, r) without and with long-distance contributions plotted in 
Figs. 3(a,c) and 3(b,d), respectively. Due to this sign alteration of Wilson coefficients in 
different SUSY models, Fig. 3 evince the variations in the magnitude of the longitudinal 
polarization asymmetries from those of the SM . This value is expected to increase in SUSY 
I and SUSY II model compared to that of the SM value due to the opposite sign Cj^C\q 
and because of the NHBs contribution in the later case. By looking at the Eq. (3.22) we 
can see that the contribution of NHBs, coming in the auxiliary functions fa and fs, is com- 
pensated by the mass of the final state lepton. Therefore the large deviation is expected 
for heavy mass of the final state lepton and large value of NHBs contribution i.e. SUSY I 
and SUSY II. Whereas, the SUSY III and SUSY SO(10) GUT models show mild effect in 
the case of muons as the final state lepton while their effects are quite distinguishable in 
the case of tauons as the final state leptons. 

Figures 4(a,c) (without long-distance contribution) and 4(b,d) (with long-distance con- 
tribution) show the dependence of normal lepton polarization asymmetries on the square 
of momentum transfer for the said decays. One can notice that the values of said asym- 
metries are quite sensitive to the large contribution of NHBs in SUSY II and SUSY III 
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models which is also clear from Eq. (3.23) whereas, it is mildly effected for the case of 
SUSY SO (10) GUT model due to small contributions of NHBs and due to the complex part 
of the Wilson coefficients. It is important to note that when there is a large contribution 




Figure 4. The dependence the probabilities of the normal lepton polarization asymmetries, Pl, 
for the decays B c — > D*£ + £~ (£ = fi, r) on q 2 without long-distance contributions (a,c) and with 
long-distance contributions (b,d), respectively, for different scenarios of MSSM and SUSY SO(10) 
GUT model. 

from the NHBs in Pjy h will changes its sign as indicated in Fig. 4 for the case of SUSY 
II. As the hadronic uncertainties are insignificant in these asymmetries so the contribution 
from different SUSY variants are quite distinguishable for the case when the final state 
leptons are muons and even more prominent when these leptons are tauons. This can 
be established from Eq. (3.23) when Pjy is proportional to the final state leptonic mass. 
Furthermore, as the normal polarization is proportional to the A which approaches to zero 
at high q 2 region and hence the normal polarization asymmetries are suppressed by A in 
this region which is depicted in Fig. 4. 

We now discuss the dependence of transverse polarization asymmetries on square of 
momentum transfer for the decays B c — > D*£ + £~ , without and with long-distance effects, 
plotted in Figs. 5(a,b) and 5(c,d) for I = fj, and r, respectively. One can see from Eq. 
(3.24), that it is proportional to the imaginary parts of the Wilson coefficients which are 
absent in the SM and in SUSY I, SUSY II and SUSY III model. But in the SUSY SO (10) 
GUT model one would expect the non-zero transverse polarization in for the B c — > D*fi + fi~ 
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Figure 5. The dependence the probabilities of 
for the decays B c -> D* s £ + £~ (£ = fj,, r) on q 2 
long-distance contributions (b,d), respectively, 
GUT model. 



the transverse lepton polarization asymmetries, Pl, 
without long-distance contributions (a,c) and with 
for different scenarios of MSSM and SUSY SO(10) 



(t + t~) decays due to complex flavor non-diagonal down-type squark mass matrix of 2nd 
and 3rd generations of order one at GUT scale, which can induce the complex couplings 
and consequently lead to complex Wilson coefficients. But these effects are very small in 
the said decay channels. These signatures of the transverse polarization asymmetries are 
depicted in Figs. 5(a,c) and Figs. 5(b,d) without and with long-distributions, respectively. 
Here we can see that its value is for small to measure experimentally both for the case of 
muons and tauons as the final state leptons. 

Apart from the above mentioned observables there is another physical observable sen- 
sitive the NP in B c — > D*£ + £~ transitions, i.e. the helicity fractions of D* vector meson 
produced in the final state. The measurement of longitudinally K* helicity fractions (//,) in 
the decay modes B — > K*£ + £~ by the BABAR Collaborations [94] put enormous interest 
in this observable. Additionally, it is also shown that the helicity fractions of final state 
meson, just like BTZ, Afb and Pl,n,t, are also very good observables to dig out the NP 
[33, 90-93]. In this respect, it is natural to study the helicity fractions for the complemen- 
tary FCNC processes such as B c — > D*£ + £~ {ft = //, r) in and beyond the SM. For this 
purpose, we have plotted the longitudinal (f£) and transverse (/t) helicity fractions of D* 
for the SM and different SUSY models in Figs. 6 and 7 for the final state leptons as muons 
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and tauons, respectively. In these graphs the values of the longitudinal {Jl) and transverse 
(/t) helicity fractions of D* are plotted against q 2 and one can clearly see that at each 
value of q 2 the sum of /l and fx is equal to one. 

Figure 6 depicts the case of muons as final state leptons, the effects of the different 
SUSY scenarios on the longitudinal (transverse) helicity fractions of D* are well distin- 
guishable throughout the q 2 region. Here one can notice that for the case of SUSY I, 
when the contributions of NHBs is neglected, the deviation from that of the SM values is 
prominent. Similarly for the SUSY II scenario, when the NHBs contributions are large and 




5 10 15 5 10 15 

<? 2 (GeV 2 ) 9 2 (GeV 2 ) 



Figure 6. The dependence the probabilities of the longitudinal and transverse helicity fractions, 
/i,Tj °f D* in B c — > D'*h + ijl~ decays on q 2 without long-distance contributions (a,c) and with 
long-distance contributions (b,d), respectively, for different scenarios of MSSM and SUSY SO(10) 
GUT model. 

tan (5 is also large, the NP effects are spotlighted. It is also clear from Fig. 6 that although 
the influence of the SUSY III and SUSY SO(10) models are mild for the case of muons 
as final state leptons but one can observe that for the case of tauons in the final state 
these effects are quite enhanced from that of the SM values (see Fig. 7). Moreover, Fig. 6 
also manifests the variations in the values of /x, {Jt) for the different SUSY variants with 
respect to that of the the SM values, which can be a good tool to put stringent constraints 
on the parameter space of different SUSY models. 

Now we turn our attention to the case, where tauns are the final state leptons, the 
helicity fractions of D* are shown in Figs. 7(a,c) and 7(b,d) without and with the long- 



- 25 - 



distance contributions, respectively. One can easily extract that similarly to the case of 





Figure 7. The dependence the probabilities of the longitudinal and transverse helicity fractions, 
/l.t, of D* in B c — > D* s t + t~ decays on q 2 without long-distance contributions (a,c) and with 
long-distance contributions (b,d), respectively, for different scenarios of MSSM and SUSY SO(fO) 
GUT model. 

muons, there is also prominent deviations in the values of the helicity fractions for all the 
SUSY models from that of the SM values. However, the effects for SUSY III and SUSY 
SO(10) models are more prominent as compare to the previous case where the muons are 
the final state leptons. These figures have also enlightened the variation in the extremum 
values of helicity fractions from the SM due to the change in the SUSY parameters. The 
deviation in extremum values are very well marked up at 12.5 GeV 2 for all the SUSY 
models, for example, the extremum value of longitudinal (transverse) helicity fraction is 
changed from its SM value 0.47(0.53) to 0.32(0.68), 0.385(0.615), 0.40(0.60) and 0.53(0.47) 
for SUSY I, SUSY II, SUSY III and SUSY SO(I0), respectively, which is suitable amount 
of deviation to measure. Hence, the measurement of the extremum values of /l and /t in 
the case of B c — > D*t + t~ can be used as a good tool in studying the NP beyond the SM 
and to distinguish among the different SUSY models. 
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5 Conclusion 



In our study on the rare B c — > D*£ + £~ decays, with I = fx, r, we have calculated branching 
ratio (BIZ), the forward-backward asymmetry Afb of the leptons, the polarization asym- 
metries Pl,n.t of final state leptons and the helicity fractions ft,T of the final state vector 
meson D* and analyzed the implications of different SUSY models on these observable for 
the said decays. The main outcomes of our analysis can be summarized as follows: 

• We have observed that the BTZs deviate sizeably from the SM value in different SUSY 
models. The study has shown that the B1Z is increased considerably for SUSY I and 
SUSY II mainly because of the change of the sign of Wilson coefficient C% . But 
for SUSY III and SUSY SO(10) the values of BTZ are mildly effected from that of 
the SM values because the small contributions of NHBs and due to the fact that 
the sign of Wilson coefficient remains the same as the SM. Hence the accurate 
measurement of the BIZs for these decays would help us to say something about the 
physics beyond the SM including SUSY models. 

• Along with the BTZ, our analysis show that Afb, especially the zero position of the 
Afb, is a fertile observable to extract the NP including SUSY. We have found that 
for SUSY I and SUSY II, Afb does not cross the zero position unlike the SM for 
B c — > D*fi + fi~ , because of the same signs of Wilson coefficient arid Cg * for 
these two scenarios of SUSY. This signature is similar to that of the observed signals 
for B —> K*fj> + fi~. Moreover, the shift in the zero positions of Afb for SUSY III 
and SUSY SO(10) models provide a promising signature of the NP which can be 
tested experimentally. Hence, the measurement of the magnitude as well as the zero 
crossing position of Afb of the considered decays can provide a stringent test for 
MSSM models and SUSY SO (10) GUT models. 

• The longitudinal, normal and transverse polarization asymmetries of leptons are cal- 
culated in different SUSY scenarios for the rare semileptonic charmed B c meson which 
can be tested in experiments with great precision. It is found that the SUSY effects 
could be measured at future experiments and will shed light on the NP signal beyond 
the SM. The transverse polarization asymmetry is the most interesting observable to 
look for the SUSY SO(10) effects where its value is non-zero in almost all q 2 region 
unlike for the cases of the SM and MSSM models. It is measurable at future experi- 
ments such as the LHC and the BTeV machines where a large number of bb pairs are 
expected to be produced. 

• We have calculated the helicity fractions Jl,t of the final state vector meson D* to 
extract the comparative study of different SUSY models and the SM. The study has 
shown that the deviation from the SM values of the helicity fractions are quite large 
with tauons in final state for the MSSM models as well as SUSY SO(10) model. It is 
also shown that there is a noticeable change due to SUSY parameter in the position 
of the extremum values of the longitudinal and transverse helicity fractions of D* 
meson for the case of tauons as a final state leptons. Hence, the helicity fraction of 
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D* meson in the decays B c — > D*£ + £ [£ = fi, r) can be a stringent test in finding 
the status of the SUSY models at the LHC. 

To sum up, the more data to be available from Tevatron and LHCb will provide a 
powerful testing ground for the SM and to put some constraints on the SUSY parameter 
space in particular the value of tan f3 and the masses of lightest chargino and lightest 
stop squark. Our comparative analysis of the SM and the different SUSY models on the 
observables for B c — > D*£ + £~ (£ = /i, r) decays can be useful for probing and its extension 
to various SUSY models the existence of the super symmetry. 
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